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Abstract. We consider a second-order equation of Duffing type. Bounds for the derivative 
of the restoring force are given which ensure the existence and uniqueness of a periodic 
solution. Furthermore, the unique periodic solution is asymptotically stable with sharp rate 
of exponential decay. In particular, for a restoring term independent of the variable f, a 
necessary and sufficient condition is obtained which guarantees the existence and uniqueness 
of a periodic solution that is stable. 

§1. Introduction. This paper is devoted to the existence, uniqueness and stability 
of periodic solutions of the Duffing-type equation 

x" + cx' + g(t,x) = h(t), (1.1) 

where g(t, x) is a T-periodic function in t and h(t) is a T-periodic function. The 
existence and multiplicity of periodic solutions of (1.1) or more general types of 
nonlinear second-order differential equations have been investigated extensively by 
many authors since C. Fabry, J. Mawhin and M.N. Nkashama initiated the study 
of the Ambrosetti-Prodi problem with periodic boundary condition [7]. However, 
the stability of periodic solutions is less extensively studied. In [18] R. Ortega 
studied (1.1) from a stability point of view and obtained an Ambrosetti-Prodi- 
type theorem under an assumption of convex nonlinearity. A.C. Lazer and P.J. 
McKenna established stability results by converting the equation (1.1) to a fixed- 
point problem [14]. Recently, more complete results concerning the stability of 
periodic solutions of (1.1) were obtained by J.M. Alonso and R. Ortega [1,2]. Under 
the condition that the derivative of the restoring force is independent of t and 
positive, they found sharp bounds that guarantee global asymptotic stability. In [1], 
optimal bounds for stability are obtained. But the above results do not cover our 
Theorem 1, since in the theorem the derivative of the restoring force may be negative 
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for some t. From their results, the key idea is to impose a condition on g(x, t) that 
can rule out the existence of additional periodic solutions that are subharmonic 
of order 2. The aim of this paper is to give conditions for existence, uniqueness, 
and local asymptotic stability. The novelty of our result is not that the friction 
constant and bounds on the derivative of the restoring force control the stability of 
the periodic solutions, but that the friction constant alone determines the rate of 
decay of the other solutions of (1.1) which are near to the unique periodic solution. 
More precisely, we will show that every solution to (1.1) that is near the unique 
periodic solution decays uniformly at the same exponential rate of c/2. Our method 
is based on linearization combined with the Floquet theory. The main results are 
the following. 

Theorem 1. Assume that g(x,t) G C 1 (R x R), and that g is T -periodic in t, such 
that 

7T 2 c 2 

(1) sap g'(t,x) = a{t) « ^ + — , 

(2) there is a T -periodic junction (3{t) G Ct such that J Q T /3(t) dt > 0, and 
g'(t, x) > (3(t) for all x G R. 

Then the differential equation (1.1) has a unique T-periodic solution, which is 
asymptotically stable. 

Here, we say that the periodic solution x of (1.1) is locally asymptotically stable 
if there exist constants C > and a > such that if x is another solution with 
||x(0) - x Q (0) II + ||x'(0) -x' (0)\\=d sufficiently small, then \\x(t) -x Q (t)\\ + \\x'(t) - 
x' Q (t)\\ < Cde~ at . The coefficient a in the exponent of this upper bound is called 
the rate of decay of x . 

Theorem 2. Assume that g(t,x) G C 1 (R x R), such that for all x G R, the 

derivative of g with respect to x is subject to the bounds 

n 2 ir 2 c 2 , . , (n + l) 2 7r 2 c 2 

-jp- + j < g x (t,x) < — — 2 + — 

for some n > 1. Then equation (1.1) has a unique T-periodic solution, which is 
stable with the rate of decay c/2. 

Theorem 3. Assume that g(t,x) = g(x) G C(R). 

n 2 ir 2 c 2 < g(x) - g(y) < (n + 1) 2 tt 2 c 2 
T 2 4" x-y ~ T 2 4 

then equation (1.1) has a unique T-periodic solution. Moreover, if h G C\ such that 
the set of critical points C — {t G [0, T] : h'(t) = 0} of h is Lebesgue-null, then the 
periodic solution is stable with rate of decay of c/2. 

Remark. The bounds in Theorem 2 are optimal for the rate of decay c/2. In the 
end of this section, we shall give an example to show that the Floquet multipliers 
associated with (1.1) may be a pair distinct real numbers, as soon as the derivative 
of the restoring force goes a little bit across the bounds given in Theorem 2. 

The following notations will be used throughout the rest of the paper. 
1. L V T : T-periodic function x G L p [0,T] with ||x||p for 1 < p < oo; 
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2. C£: T-periodic function x E C fe [0,T], k > 0, with C fc -norm; 

3. a(t) > /?(*): if a(t) > f3(t) on [0,T] and a(t) > f3(t) on some subset of 
positive measure. 

We will finish this section §1 by showing that the bounds given in Theorem 2 are 
optimal for rate of decay c/2. 
Consider the linear equation 

x" + cx' + q(t)x = 0, (1.2) 



where 

«(*) 



\ + {w + e) 2 for < t < tt, 
^ + (w - e) 2 for tt < t < TT. 



By the transformation y(t) = e ct/>2 x(t), the damping term can be eliminated and 
equation (1.2) reduced to the more familiar form of the Hill equation, 

y" + q(t)y = 0. (1.3) 

Evidently, the nontrivial solution x(t) of (1.2) has the rate of decay c/2 if and only 
if y(t) is a bounded nontrivial solution of (1.3). 

The fundamental solutions of the equation y" + w 2 y = are x\ — cos wt, X2 = 
smwt/w. Let T — 2tt, the monodromy matrix associated with (1.3) is A = A2 ■ A\, 
where w\ = (w + e) 2 and W2 = (w — e) 2 , and Ai is defined by 

( COS TTWi SmTTWi/Wi 

1 \ — Wi sin TTWi cos nwi 
and the discriminant function is given by 

( W\ W2\ . 

tr A ~ 2 cos ttw\ ■ cos TTW2 — 1 sm ttwi sm TTW2 ■ 

\W2 Wl J 

The boundary of the zone of stability of (1.3) is determined by | tr A(e)\ — 2. Ac- 
cording to [4] on page 120 (or by a straightforward computation) it can be expressed 
asymptotically by 

w = fc± £_ + o(e2) or w = fc+ i±_^_ + o(e) . 

If we choose (w,e) in the parametric resonance region, see Figure (101) in [4], 
such that (w,e) is near (n/2,0), where n is an integer, then one of the Floquet 
multipliers of (1.3) is greater than one and another is less than one. In this case, 
the trivial solution of (1.2) is still asymptotically stable, but the rate of decay is 
different from c/2. Thus if the derivative of restoring force crosses a little bit over 
the bounds given in Theorem 2, the conclusion of Theorem 2 no longer holds. Hence 
the bounds are optimal for rate of decay c/2. 

This paper is organized as follows. In §2 we recall some basic results about 
topological methods and prove a few lemmas that are crucial for the proofs of the 
main results. The proofs of Theorems 1-3 are given in §3. 
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§2. Linear periodic problems. In this section we shall recall some basic results 
about topological methods. Consider the periodic boundary value problem 

x' = f(t,x), (2.1) 

x(0) = x(T), (2.2) 

where /: [0, T] x R" — * R™ is a continuous function and T-periodic in t. In order 
to use a homotopic method to compute the degree, we assume that h: [0, T] x R™ x 
[0, 1 ] — > R" is a continuous function such that 

h{t,x, 1) = f{t,x), 

h(t,x,0) = g(x), 

where g(x) is continuous. The following continuation theorem is due to J. Mawhin 
[17]. 

Lemma 2.1. Let ft C Ct be an open bounded set such that the following conditions 
are satisfied. 

(1) There is no x £ dfl such that 

x' = h(t,x,X) VAe[0,l). 

(2) deg(.g,finR™,0) ^ 0. 

Then (2.1)-(2.2) has at least one solution. 

Let us consider the Lienard equation 

x" + f(x)x' +g(t,x) = h(t), (2.3) 

where h(t) E Ct- Evidently the periodic solution of (2.3) is equivalent to the planar 
system 

x' = y - F(x), 

y' = h(t)-g(t,x), (2.4) 

where F(x) is a primitive of f(x). A natural choice for the homotopy in applying 
Lemma 2.1 is to take 

h(t, x, y, A) = {y - F(x),\h(t) + (1 - X)h - [Xg(t, x) + (1 - A)^(a:)]}, 

where g(x) — ^ J Q T g(t,x) dt is the average of g(t,x). Since h(t,x,y,0) — (y — 
F(x), h — g{x)) = G(x), then the condition (2) in Lemma 2.1 reduces to 

dcg(G,ftnR 2 ,0) ^ 0. (2.5) 

Next we consider the system (2.1) for n = 2. We denote by x(t, x ) the initial-value 
solution of (2.1) and introduce the Poincare map P: x a — ► x(T, x a ). It is well known 
that x(t, Xo) is a T-periodic solution of system (2.1) if and only if x is a fixed point 
of P. If x is an isolated T-periodic solution of (2.1), then x is an isolated fixed 
point of P. 
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Definition. A T-periodic solution x of (2.1) will be called a nondegenerate 
T-periodic solution if the linearized equation 



does not admit a nontrivial T-periodic solution. 

Let M(t) be the fundamental matrix of (2.6) and \x\ and \ii be the eigenvalues 
of the matrix M(T). Then x(t, xo) is asymptotically stable if and only if \m\ < 1 
(i = 1, 2). Otherwise, if there is an eigenvalue of M(T) with modulus greater than 
one, then x(t, Xo) is unstable. 

Consider the homogeneous periodic equation 



where c is constant and a(t) G Lt- From now on, without further mention, we 
always suppose that the frictional constant c > 0. 

The following lemma is due to R. Ortega [18]. 

Lemma 2.2. Assume that 



Then equation (2.7) does not admit negative Floquet multipliers. 

The next lemmas are crucial to the argument for stability. 

Lemma 2.3. Assume that a{t) £ Lt satisfying (2.8) such that a(t) > 0. 

Then the moduli of the Floquet multipliers associated with equation (2.7) are less 
than one. In other words, the trivial T-periodic solution is asymptotically stable. 

Proof. We consider the following two cases. 

Case 1. If the multipliers are a pair of conjugate numbers, the conclusion of 
Lemma 2.3 follows immediately from the Jacobi-Liouvillc formula. 

Case 2. If the Floquet multipliers are real numbers, then Lemma 2.2 rules out 
negative multipliers, so it is sufficient to show that the moduli of the positive mul- 
tipliers arc less than one. If x(t) vanishes at some to, then it must vanish at T + to- 
Thus x(t) solves the boundary condition problem of the following equation: 



y' = fx(t,x)y 



(2.6) 



L a x = x" + cx' + a(t)x 



(2.7) 




(2.8) 



(e ct x')' + {a{t)e ct )x = 



x{t ) = x(T + 1 ) = 0. 



(2.9) 



Next, we consider the B.V.P. 



(e c V)' + Xe ct x = 0, 



x{to) 



x(t +T) = 0. 



(2.10) 



It is easy to verify that the n-th eigenvalue of (2.10) is 




and 
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is the first eigenfunction corresponding to the first eigenvalue 

IT 2 C 2 



We compare y(t) with the solution x(t) of the B.V.P. (2.9). 

Since a(t) <C Ai, it follows from the Sturm comparison theorem that y(t) has 
a zero in (to,T + to), a contradiction. Thus x(t) does not change sign. We may 
assume that x(t) > 0, since x'(T) = px'(0) and = Dividing (2.7) by 

x(t) and integrating by parts gives that 



r x'(t) 2 r 

/ — ny^ + clnp+ / a(t)dt = 0. 
Jo x \t) Jo 



The hypothesis of the lemma implies that lnp < 0. Hence < p < 1. 
Lemma 2.4. Assume that there is an integer n > 1 such that 

n 2 n 2 c 2 , , (n + 1) 2 tt 2 c 2 
« a(t) « y —— 4 + 



T2 4 w T 2 4 ' 

T/ien (2.7) does not admit real Floquet multipliers. 

Proof. If the conclusion does not hold, then there is a real Floquet multiplier p 
and a nontrivial solution x(t) such that x(t + T) = px(t). Since a(t) 3> Ai, it 
follows from Sturm's separation theorem that x(t) has a zero t G [0, T]. Thus 
px(to) = x(T + t ) = 0, i.e., x(t) is a solution of the B.V.P. (2.9). The assumption 
of the lemma implies that A„ a(t) <C X n +i- Therefore x(t) = 0, a contradiction. 

Lemma 2.5. Under the condition of Lemma 2.4, the rate of decay of any nontrivial 
solution of (2.9) is c/2. 

Proof. Consider the corresponding system 

X'(t)=A(t)X(t), (2.11) 

where the column vector function X{t) — (x(t), x'{t)) T and A(t) is the matrix 
function 

I 

-P(t) 



A(t) = 



Let M(t) be a fundamental matrix solution of (2.11). It is well known that M(t) 
has the form 

M{t) = P(t)e Bt (2.12) 

where P(t) and B are 2x2 matrices, P(t) = P(t + T) and B is constant. Let 
pi = e TAl and p 2 = e TX2 be the Floquet multipliers, so that Ai and A2 are the 
Floquet exponents associated with p\ and p 2 . Let x\ and xi be the eigenvectors of 
the matrix e TB . It follows from Lemma 2.2 that p\ and p 2 are a complex conjugate 
pair. Thus the eigenvectors that are associated with different eigenvalues are linearly 
independent. Therefore yi = pi{t)e Xit (for i = 1,2) form the fundamental solutions 
of equation (2.11). On the other hand, by applying the Jacobi-Liouville formula 
we have 

\ P i\ 2 =PiP2 = e-ti cdt = e- cT 
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and 

Re Ai = ReA 2 = -Rc(Ai + A 2 ) 

= ^>(w> 2 ) = -f. 

Since every solution is a linear combination of yi(t) and y 2 (*) , Pi(t) is T-periodic, 
hence it is bounded. Therefore every nonzero solution of the equation (2.11) decays 
at the same exponential rate of c/2. 

§3. Proof of main results. Now we prove our main results. 
Proof of Theorem 2. We divide the proof into three steps. 

Step 1. Existence. Without loss of generality, we may assume that g(0,t) = 0, 
for otherwise we can subtract g(0,t) from both sides of equation (1.1). Consider 
the parametrized equation 

x" + ex' + rg(t, x) + (1 - r)ax = rh(t) (3.1) 

for some a G (A„, A„ + i). We claim that there is an R > such that equation (3.1) 
has no solution on 8Br for all r G [0, 1 ]. If there is not such an R, let x n be a 
sequence such that ||ar n || — ► oo and r„ G [0,1], and denote by z n the ratio -[jf 2 ^- 
Dividing (3.1) by ||x„||, then multiplying by ip(t) G Cf. and integrating by parts, 
we have that 

f T ZnV" - ctp' + [T n g(t, x n ) + (1 - T n )a]ip dt = T f T ^ 

JO \\ x n\\ n Jo \\ x n\\ 

The condition of Theorem 2 implies that {[T n g(t, x n ) + (1 — T„)a]x„/||x n ||} is 
bounded. It is pre-compact in the weak* topology in L 1 [0, T]. Thus there is a 
subsequence such that g(t,x n )/x n — ► (3{t) and r„ —> r. Taking the limit in the 
equation (3.2), one obtains that 

/ {zip" - cz(p' + w(t)ipz} dt = (3.3) 
Jo 

where w(t) — r[3{t) + (1 — r)a. Since w(t) satisfies the condition of Lemma 2.4 
and z(t) is a T-periodic solution, it follows from Lemma 2.4 that z(t) = 0, which 
contradicts ||z(t)|| = 1. Next, by applying the homotopic invariance property, we 
have that 

deg(0 o , B R 0) - deg(0i, B R , 0) = 1. 

This completes the proof of existence. 

Step 2. Uniqueness. Let x\ and xi be two distinct T-periodic solutions of the 
equation (1.1), and x — x\ — x-i. Then x satisfies the equation 

x" + ex' + p(t)x = 0, (3.4) 

where p(t) — \g(t,X\{t)) — g{t,X\(t))\/[x\ — x^]- Since x(t) is a nontrivial T-periodic 
solution of (3.4), it will imply that x is an eigenfunction associated with a Floquet 
multiplier equal to one. Again, Lemma 2.4 rules out this possibility. Therefore 
x = 0, a contradiction. 
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Step 3. Stability. Let x(t) be the unique T-periodic solution obtained by Step 2. 
Consider the linearized equation 

v" + cv' + g'(t,x(t))v = 0. 

The condition of Theorem 2 implies that p(t) = g'(t, x(t)), verifying the hypotheses 
of Lemma 2.3, which gives the stability results about the linearized equation. 

To show that every solution of the nonlinear equation (1.1) locally decays to the 
unique T-periodic solution with rate of decay c/2, we need the following C 1 version 
of the Hartman-Grobman theorem [12]. 

Lemma 3.1. Let U be an open neighborhood of 0, and f:U C R" — ► R™ be a 

C 1 function such that /^(0):R™ — > R™ is a contraction mapping. Then f C 1 is 
conjugate equivalent to /^(0). 

Remark. The proof of the theorem depends on a C 1 version of the Hartman- 
Grobman theorem. More detailed results concerning smooth equivalence can be 
found in [6], [19], [21], [22]. In general, a C 1 hyperbolic map is only topologically 
conjugate equivalent to the linear part in a small neighborhood of the fixed point, 
thus Lemma 3.1 no longer holds for C 1 maps near the hyperbolic fixed point without 
assuming that / is contracting. 

Now we will complete the proof of Theorem 2. 

Consider the planar system associated with equation (1.1), 

x' = y — cx, 

y' = h(t)-g(t,x). (3 ' 5) 

Let Xo(t) = {xo(t),yo(t)) be the unique T-periodic solution determined by the 
initial condition Xo(0) = (xo,yo)- Then Xq corresponds to the unique fixed point 
of the Poincare mapping PX — U(T, X), where U(t, X) is the initial- value solution 
of (3.5) with U(0,X) = X. Let M(t) be the fundamental matrix solution of the 
linearization 

X' = A(t)X (3.6) 

of (3.5), where 

By the differentiability of X(t) with respect to the initial value, the Poincare map- 
ping can be expressed in terms of the initial value X by the following formula: 

PX — X = M(T)(X - X ) + o(X - X ). (3.7) 

Referring to Lemma 2.4, M(T) has a pair of conjugate eigenvalues A, A with |A| = 
e -cT/2_ xhus P(X) is a contracting mapping. According to Lemma 3.1, there 
is a C 1 diffcomorphism ip which is near enough to the identity that PX — Xq is 
conjugate equivalent to M{T). There is an invertible constant matrix C such that 

C-'M{T)C=(^ °)=£(A), 
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and we may suppose that 

±\X- X \ < \<p(X) - <p(X )\ < 2\X - X \ (3.8) 

for X — X small, since ip is near the identity. Therefore, the Liapunov exponent is 
given by 

Mx= lim -±=ln\P n X-X \ 

n—>oo nl 

= lim —\n\(poM{T) n oip- 1 (X)-ipoM(T) n o( p - 1 (X )\ 
= lim ^inl^A^oC- 1 ^- 1 ^)-^ 1 ^)]^-^ 

n—>oo nl 2 

The third equality follows from (3.8) and the fact that C is invertible. 

Remark. The above result shows that the Liapunov exponent is invariant under a 
C 1 conjugate transformation. From the proof, the conclusion is still true if ip is a 
Lipschitz mapping. 

Hence, the rate of decay of the solution to the unique T-periodic solution is c/2, 
independently of the initial value X. Theorem 1 can be proved by exactly the same 
method. 

Proof of Theorem 3. In order to prove Theorem 3, let us consider the following 
Lienard equation: 

x" + f(x)x' + g{x) = h(t). (3.9) 

The following lemma seems to be well known; however, we have been unable to 
find a published reference. Some similar results concerning sufficient conditions for 
existence were given in [9], [10]. We give a proof here, for completeness. 

Lemma 3.2. Let f,g G C(R) such that \f(x)\ > k for some k > and g(x) is 
increasing. Then (3.9) has a T-periodic solution if and only if h G ,g(R). 

Proof. The trivial necessary condition for existence can be obtained by integrating 
(1.1) and applying the mean value theorem, and the condition turns out to be 
sufficient. In fact, we consider the parametrized equation 

x" + f(x)x' + g{x) = Xh(t) + (1 - X)h = h x {t), (3.10) 

where A e [0,T]. 

Multiplying (3.10) by x' and integrating, we have that 







f \x'fdt< 


[ x'h x {t)dt 


h 


Jo 


ives that 








\Wh < 


- k \\h\\ 2 . 



(3.11) 

On the other hand, integrating (3.10) shows that there is a t G [0, T] such that 
g(x(r)) = h. Since g(x) is increasing, hence x(t) = ,9 _1 (/i) = c is unique, indepen- 
dently of A, and then a C° bound can be obtained from the following formula: 

<|c| + ± W2 . (3.12) 



<t)\ 



x{t) + x (s) ds 



It is easy now to obtain a C 2 bound for equation (3.10). 
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Let n and r-i be sufficiently large, and set 

n = {{x,y)\ M < r i, \y\ < r 2 }• 

It follows from the estimates obtained above that the equivalent planar system 
defined in Section 2, 

x' = h(t,x,y, A), 

has no solutions on dSl for A G [0, 1 ], and the computation of the degree for r\ and 
ri large enough is given by 

deg(G, n, 0) = deg( 5 (x) - h, (-n, r 2 ), 0) = 1. (3.13) 

By applying Lemma 2.1, we have that the equation (3.9) has at least one T-periodic 
solution. 

The necessary and sufficient conditions for existence of a periodic solution can 
be obtained by Lemma 3.2. In order to obtain the stability result of Theorem 3 the 
following lemma concerning the regularity of solutions of the initial-value problem 
of (1.1) is needed. Roughly speaking, under Lipschitz nonlinearity, the solution of 
the initial- value problem of (1.1) is still smooth provided that the force term h does 
not oscillate too violently. 

Lemma 3.3. Let u(t,^,rj) be the solution of the initial-value problem 

(u" + cu' + g(u) = h(t), 

\ u(o) = e, «'(o) = v ( ' J 

Assume that g is a Lipschitz function, and that h <G C\ such that the set of critical 
points C = {t G [0,T] : h'(t) = } of h is Lebesgue-null. Then, for t G [0,t], 
the partial derivatives of u and u' with respect to £, T] exist and are continuous. 
Moreover, if 

I du du \ 

<9£ dn 
du' du' ' 

V ~dj ~d^ / 



X(t) 



then 
where 



X'(t) = A(t)X(t), X(0) = Id, 
A{t)= (-g'(u(t)) -c 



Proof. The lemma says under Lipschitz nonlinearity the solution of (3.14) is still 
smooth provided that the force term h(t) does not oscillate violently. 

We were inspired by Lazer and McKenna's result concerning the regularity of 
solutions of (1.1) for the case g(u) = au + — bu~ in [15], for the general case of g(u) 
a Lipschitz function; the proof is convoluted and requires delicate analysis. We shall 
divide our proof into three steps, since it involves a great deal of real analysis. 



Step 1. 
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Claim. Let u(t) be a solution of (3.14). Let B C u([0,f]) be null. Then under the 
assumptions of Lemma 3.3, the pre-image A = u~ x (B) is null in [0,t]. 

The following result concerning measure theory is needed. 

1. Let f be differentiate on [a, b], and let A be a measurable subset of [a,b]. If 
m(f(A)) = 0, then f'(x) = a.e. x e A. 

Setting 

E n = ( x e A : ^I^hiipi > I, whenever \y - x\ < - }, 
{ \y-x\ n n J 

it is easy to check that 

oo 

E = {xe A: \f'(x)\ >0}= |J E n . 

n=l 

In order to show that E is null, it is sufficient to show that E n is null for any n. 
Now, for n fixed, by the additive property of the measure, it is enough to show that 
for any small interval I with length |/| < i, the measure m(I n E n ) = 0. Setting 
B = I n E n , since f(B) C f(A), by assumption m(f(B)) = 0. By the definition of 
a null set, there is a sequence of intervals such that 

oo oo 

U 4 =>/(£), J2 Ik<e - 

fe=l fe=l 

Setting B k = n B, then B k a E n C\I and (JfeLi s fe = B - Noting that 

f(Bk) C Ife, we have that 

oo oo 

m*(B) < ^TO*(S fe ) < ^diam(B fe ) 
fc=i fe=i 

oo oo 

< Y^ndiam(/(Sfe)) < nVm(7t) < ne. 
fc=i fe=i 

Since e is arbitrary, it follows that m(I n _E„) = 0. Hence, m(E) = 0: namely, 
for almost all i £ 4, we have f'(x) = 0. The proposition can be obtained also by 
directly applying the one-dimensional area formula. According to Theorem 3.2.5 in 
[8] on page 244 or in [16] on page 106, for a given mapping, the area of its image 
and its derivative are related by the following formula: 

2. Let f: R — > R be a Lipschitz mapping, and for any measurable set A, let N(f \ 
A,y) be the cardinal number of f^ 1 (y) D A. Then N(f\A,y)£ L X (R) 7 and 

[ \f'(x)\dx= [ N(f\A,y)dy. (3.15) 

The cardinal number N( f \ A, y ) ^ if and only if y e f(A) , so (3.15) combined 
with this assumption gives that 



f \f'(x)\dx= [ N(f\A,y)dy = 0. 

J A J f(A) 
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Thus, for almost all x € A, we have f'(x) = 0. 

We are now in a position to finish proving the claim. 

Suppose on the contrary that m(u^ 1 (Bj) > 0. It follows from above argument 
that u'(t) = a.c. t 6 A = u^ 1 (B). We may assume that A consists of accumulation 
points only, since the isolated points of A are countable, hence it is null. By Rollc's 
theorem, any accumulation point of zeros of u'(t) must be a zero of u"(t). Thus, 
u"(t) = a.c. t E A, so the Duffing's equation (3.14) reduces to 

g(u(t)) = h(t), u'(t) = a.e. t E A. (3.16) 

The above equation gives that h'(t) = a.c. t € A, because 

\g[u(t + s)] - g[u(t)]\ <L \u(t + s)-u(t)\ 
N " |*| 

where L is the Lip constant of g. This contradicts our assumption that the set of 
critical points of h is null because of m(A) > 0, which completes our claim. 

Step 2. The existence of partial derivatives. 

Let L be the Lipschitz constant of g and let {£n}i° be a sequence of numbers 
such that £ n — > £o as n — > oo. Let 

u n = u(t,£ n ,ri), u = u(t,£ ,v) 

and let 

*n = (Un - Mo)/(6i - Co)- 

Since u" + cu' n + g(u n ) = h(t) for n > 0, we have that 

|*„(t)"| <c|*„(t)'|+L|* n (t)|. (3.17) 
Moreover, we have that *„(0) = 1, *^(0) = 0, so for n > 1 and t > 0, 

*„(*) = 1+1 *„(*)' dt, K(t) = 1 + f * n (*)" dt. (3.18) 
Jo Jo 

Let M = 1 + c + L; then (3.17) gives that 

|*„(t)| + \%(t)\ < 1 + M f |*„( S )| + |<( S )| da. 

Jo 

It follows from GronwalPs inequality that 

l*n(*)| + \K(t)\ < exp(Mt). (3.19) 

From (3.17) and (3.19) it follows that the sequences and {*^(t)} are 

cquicontinuous and uniformly bounded on [0, t], so there exists a subsequence [still 
denoted by {*n(i)}] and a C 1 z(t) such that ^ n (t) — > z and *^(t) — > z' as n — > oo 
uniformly on [0,t]. Moreover, {^ n (t)} satisfies the following equation: 

*„(*)" + c*„(t)' + [gK( y'; KW)] = 0. (3.20) 

Sn - Co 
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Since g is a Lipschitz function, it follows that g'(u) exists for almost every u. Let 
E be a set such that the derivative of g does not exist; then m(E) = 0. Since 
^n(t) — > z, thus, the limits of the last term in (3.20) exist for u (t) G u ([0,t])\E. 
Let A = Uq 1 (E) be a pre-image of E. It follows from step 1 that A is null, namely, 
for almost all t G [0,t] the limit 

lim WMt)) - g(uo(t))} = j {Mmt)t 

and the Lipschitz condition gives that the above sequence is bounded. Taking 
the limit in (3.18), it follows from the Lebesgue bounded dominated convergence 
theorem that 

z(t) = l+[ z(s)ds, z'(t) = -[ cz'(s) + g'(u (s))z(s)ds. 
Jo Jo 

Thus, z'{t) is precisely the function satisfying 

z" + cz' + g'(u (t))z = 0, z(0) = l, z'(0) = 0. 

Since this determines z uniquely, the original sequences {^ n (t)} and {^(i)} must 
converge to z and z' respectively, which gives the existence of du/dt; and du'/d£. 
Similarly, we can show that du/dr] and dv! jdr\ exist. 

Step 3. Continuity of partial derivatives. 

In order to show the continuity of the partial derivatives with respect to the 
initial value, it is sufficient to show that if {?7n}i° are sequences such that 

in -> £o, Tin ->■ r?o, if u n = u(t, f„, /?„), w = u(i, £o, ??o) and if 

?y" + cy'n + g\u n {t))y n = 0, y„ = 1, ^(0) = 0, 
then y n —>■ ya, y' n — ► y uniformly on [0, f]. 

First, let us show that g'(u n {t)) — ► g'(uo(t)) in L 1 . 

Since g r (u) = lim n ^ GO n[g(u + i) — (#(«))] is the limit of a continuous function, it 
follows that g'(u) is a measurable function on X = Uo([0, i]). By Lusin's theorem, 
[8] on page 76, for every 6 > there is a bounded closed subset £1 of X = ito([0, i]) 
with m(X \ E\) < 5 such that g'(u) is continuous on E\. Again, by applying 
Lusin's theorem to the subset (X \Ei), there is a closed set E 2 of (X \E\) with 
m(X\(EiUE 2 )) < 5/2 such that g'(u) is continuous on E 2 . Evidently, E\C\E 2 = 0. 
Thus, the distance d(E\,E 2 ) > 0, which implies that the two sets do not have a 
limit point in common, so g'(u) is continuous on E\{JE 2 . By repeating this process, 
we obtain a sequence of closed subsets Fk of X such that F\ = E\ , F 2 = E\ U E 2 , 
■••7 Fk — \Jj=iEj, with the following properties: 1) g'(u) is continuous on Fk] 
2) F k C m(X \ F k ) -> 0, fc -> 00. Next, set 

G fe = { t G [ 0, t ] : u (t, Co, »») G X \ F„ }. 

We claim that m{Gk) — » as n — > 00. If not, since D Gfc+i, there is an e > 
such that m(G n ) > £o- By the monotone property of the measure, we have that 




= lim m(Gk) > £o- 

k^oo 
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But the image of A = P|fc=i° ^fc under uq is contained in P|fe=r '(X\Fk), which is 
a null subset of X. According to step 1, m(A) = 0, which leads to a contradiction. 

It follows from the theorem concerning continuous dependence on the initial value 
that u n — > Mo uniformly on [ 0, t } . Now since m{Gk) — ► as n — > oo, we may choose 
k large enough that m(Gk) < e. Now we fix fc, and since g' depends continuously 
on Fk, it follows that g'(u n (t)) — > g'(u (t)) uniformly on G% = [0,t] \ Gk, which in 
turn implies that 

lim / \g'(u n (t))-g'(u (t))\dt 



= lim / \g'(u n (t))-g'{uo(t))\dt+ [ \g' \u n {t)) - g' '{u (t))\ dt 

n ^°° JG k Jg% 

< + 2Le. 

Therefore, g'{u n (t)) — ► g'(u (tj) in the sense of L 1 . 

The rest of the proof is similar to that of Lemma 2.2 in [15]. Setting v n (t) = 
y n (t) - x n (t), then 

f v'n + ™' n + g'{u n {t))v n = [g'(u (t)) - g'(u n (t))}v' , 

\v n (0) = v' n (0). 1 ' J 

Denote the L 1 norm of g'(u n (t)) — g'(uo(t)) by e n and the C° norm of vo by fc. 
Then from (3.21), we have that 

K(t)\ + K(t)\ < e n k + M f[\v n {t)\ + \v' n (t)\] dt, 

Jo 

where M = 1 + c + L. The Gronwall inequality implies that 

\v n (t)\ + \v' n (t)\<s n ke Mt . 

This completes the proof that y n — > yo, y' n — > y' uniformly on [0,i]. 

On the same lines as the proof of Theorem 2, the stability result of Theorem 3 
follows from Lemma 3.3. 

Acknowledgment: We wish to thank the referee for his/her helpful comments 
and suggestions. 
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